Two-Dimensional Electron Gas with Cold Atoms in Non-Abelian Gauge Potentials 
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Motivated by the possibility of creating non-Abelian fields using cold atoms in optical lattices, we explore 
the richness and complexity of non-interacting two-dimensional electron gases (2DEGs) in a lattice, subjected 
to such fields. In the continuum limit, a non-Abelian system characterized by a two-component "magnetic flux" 
describes a harmonic oscillator existing in two different charge states (mimicking a particle-hole pair) where the 
coupling between the states is determined by the non-Abelian parameter, namely the difference between the two 
components of the "magnetic flux". A key feature of the non-Abelian system is a splitting of the Landau en- 
ergy levels, which broaden into bands, as the spectrum depends explicitly on the transverse momentum. These 
Landau bands result in a coarse-grained "moth", a continuum version of the generalized Hofstadter butterfly. 
Furthermore, the bands overlap, leading to effective relativistic effects. Importantly, similar features also char- 
acterize the corresponding two-dimensional lattice problem when at least one of the components of the magnetic 
flux is an irrational number. The lattice system with two competing "magnetic fluxes" penetrating the unit cell 
provides a rich environment in which to study localization phenomena. Some unique aspects of the transport 
properties of the non-Abelian system are the possibility of inducing localization by varying the quasimomen- 
tum, and the absence of localization of certain zero-energy states exhibiting a linear energy-momentum relation. 
Furthermore, non-Abelian systems provide an interesting localization scenario where the localization transition 
is accompanied by a transition from relativistic to non-relativistic theory. 

PACS numbers: 71.30.+h, 03.75.Lm, 64.60.Ak 
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I. INTRODUCTION 

Methods for creating fields that couple to neutral atoms in 
the same way that electromagnetic fields couple to charged 
particles have created the exciting possibility of studying the 
effects of a generalized magnetism using cold atoms. |lll|3,l3|l 
Using laser induced hopping, a controlled phase can be im- 
posed upon particles moving along a closed loop in an opti- 
cal lattice. The associated synthetic fields can be sufficiently 
strong to enter the regime of exotic magnetic phenomena that 
have been difficult to explore in condensed matter experi- 
ments, such as the fragmented fractal spectrum of a two di- 
mensional electron gas (2DEG) in a magnetic field, the fa- 
mous "Hofstadter butterfly" 131 ■ Such fields need not obey 
Maxwell's equations, thus providingthe possibiUty of discov- 
ering fundamentally new physics. 0!1. For example, we dis- 
cuss here the generation of non-Abelian fields, by using cold 
atoms that occupy two Zeeman states in the hyperfine ground 
level yjl ; these two states may be thought of as "colors" of the 
gauge fields, and such a system may be used to simulate lattice 
gauge theories in (2h-1) dimensions. Other potential applica- 
tions of non-Abelian fields are the creation of counterparts of 
magnetic monopoles 11'], and integer and fractional quantum 
Hall effects. H 

In this paper, we adopt the 2DEG as a motif for the study 
of cold atom systems. The subject of 2DEGs in a magnetic 
field is a textbook topic ||6'1, as the problem maps to a one- 
dimensional harmonic oscillator. The discrete energy lev- 
els of the oscillator are the Landau levels that describe free 
particle energies in terms of the quantized units hojc, where 
ujc = eB/mc is the cyclotron frequency of the corresponding 
classical motion. Each level is highly degenerate, reflecting 
the fact that a classical electron spirals about a line parallel to 



the magnetic field, with an arbitrary center in the transverse 
plane. The degree of degeneracy is equal to V^ /2tt5'^ where 
5 ~ ^Jhc/eB is the magnetic length and L? is the area of the 
system. 

Beginning with the celebrated work of Onsager Q, 
Harper jS] and then Hofstadter 121, the subject of 2DEGs in 
a crystalline lattice in a magnetic field has fascinated physi- 
cists as well as mathematicians. In the presence of a lattice, 
each Landau energy level splits into Q bands, where the ra- 
tional number P/Q is the magnetic flux through the unit cell 
in units of the magnetic flux quantum (fluxoid). The heart of 
the problem is the two competing periodicities related to the 
ratio of the reciprocal of the cyclotron frequency and the pe- 
riod of the motion of the electron in the periodic lattice. Two 
key aspects that have been explored extensively are the exotic 
multifractal spectrum (Hofstadter butterfly), and the metal- 
insulator transition obtained by tuning the ratio of the tunnel- 
ing along the two directions of the lattice I8f|. Recent studies 
have shown that these properties can be demonstrated using 
ultracold atoms in an artificial magnetic field 12, [i3l • This 
paper revisits the metal-insulator transition when the 2DEG is 
subjected to a non-Abelian gauge field which is a natural gen- 
eralization of the uniform magnetic field. Such fields yield a 
much richer spectral and transport landscape than is encoun- 
tered in the Abelian case. 

The generic experimental setup for producing non-Abelian 
U{2) gauge fields that we consider here, consists of a two- 
dimensional optical lattice populated with cold atoms that oc- 
cupy two hyperfine states 13, [211 ■ Such systems exhibit three 
competing length scales , associated with two distinct "mag- 
netic fluxes" (denoted by ai and a2) that penetrate the unit 
cell. Our aim is to describe some of the generic properties 
of such systems. Although our main focus is on optical lat- 



tices, we first discuss the corresponding continuum problem, 
where the the infinite degeneracy of the Landau levels is lifted 
by non-Abelian interactions. The continuum problem mirrors 
some of the features subsequently encountered in the lattice 
system. 

In the discussion of the metal-insulator transition in the lat- 
tice, we focus on the ground state as well as the states at the 
band center. These two cases are respectively relevant for ex- 
perimental systems involving Bose condensates and fermionic 
gases near half-filling. Some of these results were described 
in an earlier paper II12I1 . In addition to a detailed analysis, here 
we describe new results, such as the simulation of relativistic 
phenomena using cold atoms in non-Abelian fields. By tun- 
ing lattice anisotropy, we can implement relativistic as well 
as non-relativistic dynamics, with a particular focus on the ef- 
fects of disorder. Simulation and detection of Dirac fermions 
using single-component cold atoms in a hexagonal lattice was 
recently proposed lllSll . The systems we propose here provide 
the possibility of observing relativistic particles and also of 
studying their localization properties. We show that the non- 
Abelian systems provide an experimental realization of the 
defiance of localization by disordered relativistic fermions, a 
topic that has been the subject of extensive study [14]. 

In Section II, we introduce non-Abelian gauge fields and 
the corresponding effective "magnetic fields." Section III ex- 
amines the continuum limit of a single particle in a non- 
Abelian gauge field. In Section IV, we discuss lattice systems 
subjected to these fields, and describe methods of calculation. 
In Section V, we study various spectral characteristics of the 
non-Abelian lattice systems. There, following long establised 
practice for studying metal-insulator transition in Abelian sys- 
tems, we fix ai = {Vb — l)/2, the golden mean, which we 
denote as 7. The irrationality of ai ensures the existence of a 
localization transition.lJSl Ulll For a2, we consider a selected 
set of both rational and irrational values. Sections VI and VII 
discuss localization properties of the states at the band cen- 
ter {E = 0) and at the band edge. The localization of the 
E = states brings out some of the most important fea- 
tures of the non-Abelian cases, including the dependence of 
the transition upon a conserved momentum. Furthermore, a 
unique aspect of the the non-Abelian system, namely the de- 
fiance of localization of the E = states, emerges when the 
energy-momentum relation mimics the behavior of relativistic 
particles. Section VIII describes the experimental realization 
of the metal-insulator transition in cold atom lattices. 



II. NON-ABELIAN GAUGE FIELDS 

Effective non-Abelian vector potentials arise naturally in 
systems where the atoms have N degenerate internal states. 
The most general vector potential couples the states, and thus 
gives rise to a U{N) gauge symmetry. We here consider the 
case where A^ = 2. In our treatment of the non-Abelian case, 
we follow the convention of an earlier study 1311, adopting its 
form of vector potential, 
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The ai determine the "magnetic fluxes" of the lattice with lat- 
tice constant a. 

Equation ^ is in the Landau gauge: A{x,y) = 
{Ax, Ay{x), 0), where A^ is a constant and Ay depends only 
on X. 

We rewrite the vector potential in terms of Pauli matrices 
(Ti, separating the Abelian and the non-Abelian parts of the 
gauge field. 



lie r TT X X 

A = — --(/-f7a.)x + 27r-a/y + A-cr^y 
2 a a 



, (2) 



where we have defined quantities a = (cvi + q;2)/2 and A = 
TT{ai — 0^2). Here aa, a = x,y,z denotes Pauli matrices. 
The parameter A characterizes the non-Abelian feature of the 
system. 

For non-AbeUan fields, the effective "magnetic field" is 
given by. 



B = V X A- —A X A. 
he 
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The origin of the extra term A x A can be traced to the com- 
mutator for the generalized velocity operator {p — -A)/M, 
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For the vector potential in Eq. ([T]), this gives 
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where Bq = 27ra(^). Thus, a = Boa^ /{2TThc/e) de- 
scribes the Abelian flux quanta penetrating per unit cell of 
the lattice. The non-Abelian gauge potential generates a non- 
uniform magnetic field, as B depends expUcitly on the spatial 
coordinate x when A 7^ 0. 



m. CONTINUUM LIMIT OF THE NON-ABELIAN 
SYSTEM 

We now consider the continuum limit of the non-Abelian 
problem. Although the A of the Eq. (|2]i is ill-defined in 
the continuum limit a ^ 0, the study is useful in illustrat- 
ing some key aspects of the non-Abelian systems. In general, 
continuum problems can also be experimentally realized, as 
inRef. H. 

It can be shown, after some algebra, that the two- 
component continuum Hamiltonian He = (p — -A)^/(2A/) 
resulting from the vector potential in Eq. (yj is gauge- 
equivalent to 



fr _ J_f{p.+f^? + V{x) CA{x-xo)' ,,,. 



up to a ky dependent term. The transverse momentum ky 
is a conserved quantity as the Hamihonian He with A given 
by Eq [2] is cyclic in y. Here, /? = ||, C = 8n^a{P/a)^ 

M^oj^xy2 with cj = ''-MVW 

hkyP aia2_ 



y{x) = 

xq = 27ra 



al)/2 and 

This particular form of the Hamilto- 
nian provides a new, illuminating picture of the non-Abelian 
problem; the particle behaves as a two-component harmonic 
oscillator existing in a positive as well as a negative charge 
state. The physics of this system is that of a particle-hole pair, 
with the non-AbeUan parameter A , governing the coupling 
between states. 

The spectrum of He is obtained by numerical diagonaliza- 
tion in a basis of harmonic oscillator wave functions with fre- 
quency ui. Figure ([T]) shows the six lowest energy levels. For 
fixed ky, at A = 0, each Landau level is two-fold degenerate. 
For A 7^ the degeneracy is lifted and the eigenstates become 
entangled states of a particle-hole pair 
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FIG. 1: Six lowest energy levels of the continuum non-Abelian 
Hamiltonian in Eq. (O, /3 = 7r/2, M = l,ky — 1 and 27rai = 1. 
The levels are equally spaced in (and only in) the Abelian case, 
A = 0, where the system reduces to two decoupled harmonic os- 
cillators. 
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FIG. 2: (color online) The ground state splitting obtained numeri- 
cally (see Fig. [Til (solid black), compared with that obtained by per- 
turbation theory (dashed red) about the Abelian point A = 0. 



Figure [3] shows the variation of the energies with ky, ob- 
tained numerically. In highly non-Abelian cases, the energies 
bear no relation to their Abelian values. Close to the Abelian 
limit (bottom) the energy levels are simply split around the 
Abelian energies. The energies oscillate with ky, resulting in 
actual and avoided crossings (i.e., the Landau bands overlap). 
In the vicinity of the crossings, the bands exhibit a linear dis- 
persion relation. As we shall discuss, these features reappear 
in the corresponding problem of the non-Abelian gauge field 
on an optical lattice. 

Figure m summarizes the effects of the non-Abelian gauge 
potential on the lowest Landau level of the coiTesponding 
Abelian problem. The figure describes the continuum limit of 
the Hofstadter "moth" f^, which is the generalization of the 
Hofstadter "butterfly" as the underlying gauge field becomes 
non-Abelian. This coarse-grained "moth" illustrates the sym- 
metry breaking feature of the non-Abelian system as it lifts 
the degeneracy of the corresponding Abelian problem. 



IV. TWO-DIMENSIONAL LATTICE IN NON-ABELIAN 
GAUGE FIELDS 



As Fig. [T] shows, the energy levels are equally spaced only 
for the Abelian case, ai ~ a2- We can explicitly understand 
the splitting of each Landau level via degenerate perturba- 
tion theory, with A as a small parameter and using degenerate 
eigenstates. 
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where ipn{x) are the eigenstates of the corresponding har- 
monic oscillator Figure [2] compares the perturbative splitting 
of the lowest Landau level with the numerical result. 



Our starting point is a tight binding model (TBM) of a par- 
ticle moving on a two-dimensional rectangular lattice {x,y), 
with lattice constants (a, b) and nearest-neighbor hopping 
characterized by the tunneling amplitudes (J, J A). When a 
weak external vector potential, A{x, y) ~ {Ax,Ay, 0), is ap- 
plied to the system, the Hamiltonian, 



H = -J 



cos ( {px - -Ar,)- ) + Acos ( {py Ay)- 



where pis the momentum operator Alternatively, the Hamilo- 
nian of a 2DEG on a lattice in the presence of a magnetic fied 
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FIG. 4: (color online) The continuum version of the Hofstadter 
"moth" |y]- This plot shows the energies as a function of qi and 
a2 for a range of ky ; the color scale indicates the range of ky. Along 
the line ai = 02 , the Abelian "backbone" of the moth, there is no 
ky dependence. 



is the magnetic flux penetrating the unit cell in units of mag- 
netic flux quantum, $0 = ch/e. 

We denote the eigenfunction (projected onto the x, y basis), 
corresponding to the eigenvalue equation iJ|vl> >= E\ip > 
as ^(.T, 1/). With the transverse wave number of the plane 
wave as ky = ky/a, the wave function can be written as: 
'i'lma, no) = e^'^'^'^^^gm with x = ma and y = na. 

Subsituting the vector potential defined in the Eq. |2] the 

two-component vector g™ := | " '" | can be shown to result 
in the following equations. 
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FIG. 3: (color online) The three lowest energy levels of the non- 
Abelian system with ai = 1.32 and a^ = 0.253 (top), 02 — 1.11 
(bottom). The dashed line shows the Abelian levels corresponding to 



where 



\/|(a?+ai). 



can be written as. 



h = -J2j, 






H.C. 



(6) 
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where c^ is the usual fermion operator at site i. The Jy is the 
nearest-neighbor anisotropic hopping with values J and JA 
along the x and the y-direction. 

The phase factor 6'y = —Oji defined on a link < i, j > is 
identified as (27re/c/i) § A.dl, where A is the vector poten- 
tial, and 
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Um = 2Acos(27raiTO — 27rfcj,), 
Vm = 2Acos(27ra2?n — 27rfcj,). 

For ai = a2 (mod 1), we recover the Abelian limit de- 
scribed by the Harper equation [Sfl 

9m+i + g-in-i + 2Acos(27ram - 2nky)gjn = E g^n, (8) 

For iiTational values of the flux a, the system exhibits a metal- 
insulator transition at A = 1. 

The approach to the irrational values of a; is studied by 
considering a sequence of periodic systems obtained by ratio- 
nal approximants a^ = Pi/qi- This corresponds to truncat- 
ing the continued fractional expansion of ai and 02. The re- 
sulting periodic system will have period Q, the least common 
multiple of qi and (72- 

The 2(5-dimensional system can be cast in the form of two 
Q-dimensional eigenvalue problems: 
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It is instructive to compare the eigenvalue formulation to 
the transfer matrix approach discussed in earlier studies y]]. 
1 \The TBM equation can be written as a transfer matrix equa- 
tion. 
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Here Ea and Eb denote the two sets of eigenvalues of 
two uncoupled systems. The allowed eigenenergies of the 
system are the union of these two sets. 

In the above two eigenvalue equations, we have used 
Bloch condition. 



The allowed energies are those for which the product of 

Q successive matrices B{m) has an eigenvalue on the unit 

circle, e^'^^^ and e^"^^. Alternatively, the 4-dimensional 

\Vi / transfer matrix equation can be reduced to two independent 

2-dimensional transfer matrices as 
the 
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where the 2x2 matrices Ta and Tg are given by, 

TA{m) = 
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FIG. 5: (color online) The two possible antiferrimagnetic states des- 
ignated as A and B. 

An important consequence of this decoupling of the 2Q- 
dimensional problem into two Q-dimensional problems is that 
the eigenstates of the system are of an "antiferrimagnetic" 
type, as shown schematically in Fig.|5] We will refer to them 
as of the A and B-type. The corresponding states denoted as 
XA and xb are in general non-degenerate. 



This decoupling of the 4-dimensional transfer matrix prob- 
lem into two 2-dimensional transfer matrices is equivalent to 
the decoupling discussed earlier for the eigenvalue problem, 
which in turn implies the possibility of "antiferrimagnetic" 
type states as shown in Fig.|5] An important consequence of 
this type of state is that (out of four), only two of the eigenval- 
ues of the 4-dimensional transfer matrix have to be on the unit 
circle. In other words, in contrast to the statement made in 
an earlier paper yl], the allowed energies include states where 
two of the four eigenvalues of the transfer matrix are not on 
the unit circle III2I1 . 

The existence of "antiferrimagnetic" states and the relation- 
ship between the direct diagonalization method and the trans- 
fer matrix approach can be illustrated by considering a sim- 
ple non-Abelian system, namely the one with ai = 1/2 and 
a2 = which can be treated analytically. 

Diagonalization of two independent 2x2 matrices Eq. ^ 
gives 

EA{kA,ky) = 2 [A cos 27rfcj, ± cos A:^] , 

EB{kB,ky) — ±2 J A'^cos'^ 2TTky +cos'^ kB, 



and the corresponding eigenvectors 



XA = 




XB = 
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In general, the A and the B-states are non-degenerate. As is 
explicit in this example, the magnitude of the two components 
of the vectors are in general unequal, and hence the solutions 
correspond to antiferrimagnetic states. However, at ky = tt/2, 
the Ea ~ Eb and the two degenerate states are antiferrimag- 
netic. 

The spectrum can also be obtained by iterating the transfer 
matrix problem where the energies can be written in terms of 
the eigenvalues of the transfer matrix, e^"^-^ , e^*"^^ . Com- 
parison of the spectrum obtained by these two methods shows 
that (t)A — kAl'2. and (^b = kB/'i- 

The results of this paper were obtained using the direct di- 
agonalzation method. We shall henceforth set ai = 7, the 
golden mean, and explore various complexities of the prob- 
lem for different values of 02. 
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V. NON- ABELIAN SPECTRUM FOR IRRATIONAL a 1 

The energy spectrum of the system is the union over ky of 
the individual energy spectra of the A and the B types of the 
tight binding equations. In the Abelian case, ai = a2, and 
when ai is rational, equal to p/q, the spectrum consists of q 
bands which are usually separated by gaps. As ky varies, the 
bands shift and their width may change, but they do not over- 
lap, except at the band edges. For irrational a, the spectrum is 
independent of ky . 

A striking aspect of the non- Abelian problem is the overlap- 
ping of the bands, as illustrated in figures|6]and |2l These fig- 
ures depict two different classes of typical non- Abelian spec- 
tra with rational a2: In Fig. |6l type A and B solutions result in 
a non-degenerate spectrum; and Fig. |2]shows the case where 
A and B solutions are degenerate. Below we discuss various 
spectral characteristics of the system. 

For rational a2 = P2/'?2, the spectrum is a periodic func- 
tion of ky. This is due to the fact that for irrational ai, the 
set {Urn} is ergodic in in, while the set {V„i} is periodic in 
■m for rational a2- We list below some of the characteristic 
properties of the spectrum: 

(1) EA,Bi2nky) = -EA,Bi2Trky + 2^2), 

(2) EA{2TTky) = -EB{2nky + ^i^a^), 

(3) EA,B{2nky) = EA.B{2nky + 477^2), 

(4) E{2Trky) = E{2Trky + 47ra2). 

For certain values of a2, A and B-type of states are de- 
generate. This happens when the two sets Vzmiky) and 
V2m+i{—ky) coincide. This degeneracy occurs when (a) a2 
is an irrational number and (b) 02 = p/q with q-odd as 
Vm+qn{ky) = Vq-rn+q7i{-ky) (whcrc n is an integer and 



FIG. 6: (color online) Energy spectrum viewed as a function of ky 
for Qi = 89/144 and a2 — \ with A = 1 for a range of kx values. 
The red and blue correspond to Ea and Eb respectively. The grey 
dots show the corresponding Abelian case with oi = 02 = 7- 
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FIG. 7: (color online) Energy spectrum viewed as a function of ky 
for Qi = 89/144 and 02 = | with A = 0.5. The grey dots show a 
corresponding Abelian case with ai — C12 — k- 



m < q). For even q, the A and B states are in general non- 
degenerate. This distinction between the odd and the even 
cases leads to significant differences between these two cases. 



VI. THE LOCALIZATION TRANSITION 

The metal-insulator transition llllll in 2DEGs in the pres- 
ence of a magnetic (Abelian) field is a paradigm for the Ander- 
son localization transition. We now discuss the corresponding 
localization transition that exists in the non-Abelian systems. 
In contrast to the Abelian case, where all states localize at 
the same value of the tunneling anisotropy, localization in the 
non-Abelian case varies throughout the spectrum. 

In this section, we will discuss the localization properties 
of the E ~ {] state a study is relevant for fermionic atoms 
near half-filling. As shown below, for 02 = ^ as well as for 
a2 = \, the onset to a localization transition can be inferred 
from the well-known localization characteristics of the Harper 
equation. 



Al: Localization Boundary for a2 = 1/2 
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FIG. 8: (color online) The eigenvalues of the effective TBM 
(Eq. dlOb ) describing E = G states as A varies. The line (red on- 
line) shows the values of A where _E = is an eigenvalue of the 
TBM(Eq.[TOj 
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FIG. 9: (color online) With 02 = 1/2 shaded regime (blue) shows 
the extended phase while shaded red shows the localized phase of 
£■ = state in A — fc^ plane. 



where e = —2 and A = 2A^ cos(27rfcy). 

For _E = 0, the uncoupled 0-equation maps to an £' = — 2 
Harper-like equation (Eq. (fTOli). where the on-site quasiperi- 
odic potential is a sinusoidal function of ky. The eigenstates 
of this system localize at A = 1, providing an explicit thresh- 
old for localization of the E ~ Q state of the non-Abelian 
system provided e = —2 is the eigenvalue of Eq. (fTol i. 

As shown in Fig.|8] e = —2 is an eigenvalue of the system 
provided A = Ai < 0.48 or A = A2 > 1.83. These critical 
values determine the boundary curves for the localization 
of the E ~ state: in the Harper equation, all states are 
extended for values of A < 1. These two localization 
boundaries are exhibited in Fig. ?? 

A2: Localization Boundary for a2 = 1/4 

For a2 ^ \, the uncoupled 6'-equations for the A and the 
B-sectors of the TBM for S = reduce to 



qA 
'^m+2 



^m-2 
^m-2 



2i\A cos(27raiTO — 2TTky)9^^ 
- 2i\B sin(27rQi77i — 2TTky)9^ 



= 0, 
= 0, 



where 



For a2 = ^, the coupled TBM equations (Eq. [TOl ) for 
E = reduce to 

0m+2 + Om-2 + 2{-l)"'Xcos{2TTaim-2Trky)0,n = eOm 

Vm+i + Vm-i + 2A coa{2iTaiin — 2iTky)dm = 



A^ = 2A^cos(27rA:y), 
As = 2A^sin(27rfc.y). 

The above two equations correspond to A and B-type states 
with E — 0, respectively. Here ^jfj''^ = i"'9^'^ . In a manner 



analogous to the corresponding Hermitian problem, the sys- 
tem exhibits self-duality at A^ ^ = 1 and this self-dual point 



(mod 



describes the onset of localization lITsll . For ky 
i), both types of solutions localize simultaneously. However, 
at other values of the transverse momentum, only one of the 
states is localized. This is an example of two degenerate states 
with different transport properties: depending upon (A, ky), 
type A states may be extended (localized) while type B states 
will be localized (extended). This localization boundary in 
A — ky space is shown for types A and B in Fig.fTol 
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FIG. 10: (color online) For a2 — 1/4, shaded regimes show ex- 
tended phase for two degenerate E = states belonging to the A 
(green) and B (blue) sectors in A — fcy plane. 

The existence of conducting states for all values of A is 
one of the most intriguing characteristics of the non-Abelian 
system. Below we show that these states defying localization 
describe relativistic particles. 

B: Relativistic Dispersion and Defiance of LocaUzation 

Figure [TT| shows the energy-momentum relation for a2 = 
1/2 near E ~ Q,ky ~ 1/4. Although the level structure 
is complicated, near ky = 1/4, the energy bands exhibit the 
linear dispersion characteristic of the one-dimensional rela- 
tivistic particles. Thus, the non-Abelian system with A and B 
type states, provides an interesting manifestation of the posi- 
tive and the negative energy states of a one-dimensional rela- 
tivistic particle. 

An important characteristic of the states that reside at the 
crossings is that they defy localization. It should be noted that 
a crossing at E ^ exists irrespective of the value of A. In 
other words, we have a relativistic theory for all values of A as 
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FIG. II: (color online) Blowup of the E = and ky — 1/4 neigh- 



i and A 
blue respectively correspond to Ea and Et- 



borhood for ai = 89/144 and a2 — 2 



10. The red and 



shown in the Figures |6] and [TT] Such states remain extended 
irrespective of the quasiperiodic disorder in the system as the 
linear dispersion exists for the full range of A values. 



For a2 



we obtain an effective relativistic theory for 



zero-energy states near ky = for type-A and near ky = 
tt/2 for type-B states. These states remain conducting for all 
values of A. 

We would like to note that in the Abelian system, a linear 
energy-momentum relation resulting in a Dirac cone occurs 
for rational values of a in the neighborhood of some spe- 
cial values of k^^ky near E = Q. However, for irrational 
a, the spectrum is independent of ky and the Dirac cone dis- 
appears. Therefore, in the Abelian case, effective relativistic 
theory bears no relationship to the transport properties as the 
states are always extended for rational a. 

C: LocaUzation Transition and Loss of Relativistic 
Dispersion 

Our detailed investigation for various values of 02 shows 
that the presence of conducting states for all values of A is not 
a generic property of the system. In particular, for cases where 
the type-A and type-B states are always degenerate, all states 
are found to localize. Interestingly, the transition to localiza- 
tion is accompanied by a loss of the relativistic character of 
the energy momentum relation. 

For example, for a2 ~ p/q where q is odd, as well as for ir- 
rational 0:2, the crossings characterizing certain E ~ Q states 
disappear beyond a certain critical value of A. Interestingly, 



this threshold for the disappearance of the crossing is always 
found to coincide with the onset to localization of that state. 
Figures [12] and [13] illustrate this for irrational 012 as the disap- 
pearance of band crossings is accompanied by the broadening 
and flattening of the Bragg peaks. 
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(color online) Fourier transfonn of the wave function for 

3 



FIG, 

E = state with 02 = 7'' with A = 0.75 (sharp fringes) and A 
1.25 (smeared out fringes). The label k on the x-axis corresponds to 
the momentum 2Trk/L, where L is the size of the lattice. 



FIG. 12: (color online) Top and bottom panel respectively show the 
spectrum for 02 — 7'^ with A = 0.75, 1.25 which respectively cor- 
respond to extended and localized E = states. 



VII. LOCALIZATION TRANSITION OF BOSE-EINSTEIN 
CONDENSATES 



The natural locus for BEC in ultracold atoms in optical lat- 
tices is the band edge. We now explore the spectral and trans- 
port properties of the states at band edges, namely the mini- 
mum energy states as ky varies. 

In contrast to the preceding analytical treatment of the band 
centers, we have investigated localization properties of the 
states at the band edge with numerical methods. 

As seen from Figs |6] and [141 the energy spectrum for 
a2 — -^ shows the existence of a linear dispersion relation 
near the band crossings. As the lattice anisotropy A varies, 
we see a transition from relativistic to non-relativistic behav- 
ior near A « 2.5; this transition is accompanied by the loss of 
the wave function's spinor character, causing an effective spin 
polarization. 

The robustness of the hnear dispersion in non-Abelian sys- 
tems is shown for various values of a2 in Figs. [15] and [16] 
It appears that it is only in the even-g cases that the nature of 
the dispersion changes as A varies. Similarly, for a2 = 7'^ 
(an odd harmonic of 7, as 7"^ = 2 — 87), linear dispersion 
at fcj, = and at ky = 1/4 occurs for all values of A, while 
for a2 = 7^ (an even harmonic of 7, as 7"^ = 27 — 1), a 
relativistic energy-momentum relation is seen for small and 
large values of A as illustrated in Fig. [16] In other words, a 
"transition" from relativistic to non-relativistic behavior can 
be induced by varying A for some values of a2- 
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FIG. 14: (color online) Minimum energy as a function of ky for 02 — 
^ with A = 0.5, 1.5,2.5,3.5 (top-bottom) illustrating the change 
from linear to quadratic dispersion near ky = ±1/4. 
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Another point to be noted is that the ground state of the 
system may have nonzero momentum: for even q, the global 
energy minimum occurs at ky = 0, while for odd q, it occurs 
at ky = a2/2. 

Figures [TTl illustrate the localization transition of the mini- 
mum energy states. Extended states in these figures are char- 
acterized by sharp Bragg peaks in the momentum distribution, 
and the localization transition is signaled by the broadening of 
these peaks. As we increase the parameter A, the ky = ±1/4 
states localize before the ky = state. Our detailed inves- 
tigation shows that ky = is the last state to localize as A 
is varied for all values of 02. This is contrary to familiar ex- 
perience, in which localization begins at the band edge. The 
localization for the minimum energy state is insensitive to the 
energy-momentum relation, in contrast to the _E = states. 

For irrational 012, we expect the localization threshold to 
be lowered. Our numerical results show that the minimum 
energy states begin to localize at a relatively small value of 
A « 0.15. As discussed earlier, E = states resist localiza- 
tion due to their linear dispersion but eventually localize. Our 
numerical studies show that localization is complete at A = 1, 
as in the Abelian case. 



FIG. 15: (color online) Minimum energy as a function of ky for 
02 = 1/4 (cross), 02 = 1/5 (dots) with A = 1.5,2.5,3.5,4.5 
(top-bottom) illustrating the change from the linear to quadratic dis- 
persion near ky = ±1/4. 



VIII. EXPERIMENTAL OBSERVATION OF 
METAL-INSULATOR TRANSITION 



An experimental setup for generating artificial Abelian and 
non- Abelian fields consists of f2|, 9] a two-dimensional opti- 
cal lattice populated with cold atoms that occupy two hyper- 
fine states. The lattice laser polarization is adjusted to confine 
these states to alternating columns. The non-Abelian scheme 
requires atoms with two pairs of hyperfine levels: \gi), |ei), 
I32), |e2) as shown in Fig.fTSJ 

The typical kinetic energy tunneling along the y-direction 
is suppressed by accelerating the system or applying an in- 
homogeneous electric field in that direction such that the lat- 
tice is tilted. Tunneling is instead accomplished with two sets 
of laser-driven Raman transitions with space-dependent Rabi 
couplings iljc'^'^^y where j ~ 1,2. The wave numbers qj gen- 
erate an effective magnetic flux where qj = (27raj)/a, where 
A = 2a is the wavelength of the laser light. In an optical 
lattice with a finite number of sites, the two components of 
the "magnetic flux" ( ai, 02 ) can be adjusted, in a controlled 
manner, to a sequence of rational approximants to the golden 
mean by tuning the qj . We direct readers to Refs. 0[2ll regard- 
ing various details for generating these artificial gauge fields. 
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FIG. 16: (color online) Variation of minimum energy with ky for 
a2 = 7^ (large crosses) and 02 = 7^ (smaller crosses) with A — 
1.5, 2.5, 3.5, 4.5, 5.5, 6.5 (top-bottom). 



We now describe the experimental feasibility of tuning A 
to induce metal-insulator transitions by adjusting the lattice 
beam intensity Vq. For simplicity, we will initially discuss the 
Abelian case. Let us first consider the laser assisted coupling 
Jy as a function of Vo/Efi, where Er = 2T:^fP/M\'^ is the 
photon recoil energy. The tunneling is defined as the matrix 
element of the Rabi coupling (il) between Wannier functions 
w, evaluated at the two adjacent lattice sites: 



w{x ~ Xi) — Q.cxp{iqx)w{x — Xi-i)d^x, (10) 



where q ~ (27ra)/a. The Wannier functions for V{x) = 
Vo sin^(27rx/A)] have been computed 131; Jy decreases 
monotonically with Vo/Efj. This basic behavior can be 
demonstrated analytically by assuming a deep lattice approx- 
imated by a harmonic oscillator potential and taking the Wan- 
nier functions to have the corresponding Gaussian form. The 



FIG. 17: (color online) Fourier transform of the wave function for 
A — .4, .5, .7 (top-bottom). Each caption shows ky = .25(red) and 
ky = O(blue) for the minimum energy state for Q2 — 1/2. The 
X-axis index k corresponds to the Bloch vector 2-Kk/L. There are 
peaks corresponding to the irrational values of oi, occuring at k val- 
ues equal to half of the Fibonacci numbers (as anti-ferromagnetic na- 
ture effectively doubles the size of the unit cell). Additional satellite 
peaks characterize the non- Abelian ferature of the system. 



Gaussian approximation yields 



hVL 



■cxp[ 



TH^^Wl^l-'l-^WSl 



(11) 



The kinetic energy coupling in the a--direction J^ also de- 
creases monotonically with Vq/Eh for sufficiently large val- 
ues ifisll as described by 



J,^ « l.-i^lER 



Er 



1.051 



exp[-2.l2l^Vo/En]. (12) 



The ratio of A = J.,/Jy = {Er/HQ) /{Vo/Er, a) is 
shown in Fig.[T9]for a characteristic range of Vq/Er^ with the 
scale set by the factor KI/Er. In order to generate a useful 
range of A values (e.^.,0 < A < 2), the parameter Ml/ Er 
must be set to order unity. 
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FIG. 18: (color online) Schematic diagram illustrating the non- 
Abelian t/(2) set up. The ground states are hollow and the excited 
states are filled. Red and blue represent the two "colors" of the C/(2) 
group. 



8r ~ 100 GHz. We require that the Rabi coupling of the 
Raman transition J7, the detuning A, and the lattice trapping 
frequency v^ = VlE^jV^/fi, have well separated magnitudes 
such that 17 ^ A ^ i/j;, to ensure that only the lowest band 
of the lattice is occupied and no other excitations occur. Typ- 
ical values of v^ are on the order of tens of kilohertz. The 
Raman transition is stimulated by two lasers with Rabi cou- 



>(i 



>(i) 



plings rig and ile and intensities Ig and /g with a large 
detuning 8r such that the effective Rabi coupling magnitude 
is O = J7g -^riji ' /25r. The Rabi coupling Vl can be written as 
a product of atomic factors and laser tuning parameters. 



^ = 






where F is the natural decay rate of the 6^^3/2° state and I gat 
is the saturation intensity of the D2 line (See Ref. [17]). The 
ratio ^ = le/Ig must be less than 0.17 or greater than 5.8 to 
satisfy the J7 <C A condition. The separation of scales be- 
tween the Rabi coupling Q, of the Raman transition and lattice 
trapping frequency v^ necessary to generate the magnetic field 
in the above scheme (i.e., Vi ^ v^.) is sufficient to generate a 
reasonable range of A values. 

In the non-Abelian case, there are generally two possible 
values of A corresponding to f^i and 0,2, one for each "color". 
By adjusting VI2/O1 = f{Vo/ER,ai)/f{Vo/ER,a2), we 
obtain a single A in correspondence with the theoretical stud- 
ies described here. 




FIG. 19: (color online) A as the depth of the 2D optical lattice is 
tuned with a = 5. The factor Er/MI Ri 4.2 with the following laser 



parameters: Ig 
3.2 kHz. 



1 mW/cm^ 



X 



11, 5r = 100 GHz, En/h 



We now argue that it is possible to achieve this with rea- 
sonable experimental settings. We consider the case of ^^Rb, 
where the \g) and |e) states are taken to be the hyperfine lev- 
els of the 5^5*1/2 level lUTIl and the Raman level is S^Pg"^- 
The parameter MI/Er can be fixed near unity if the Raman 
laser beams have intensities on the order of 1 mW/cm^ with 
a; « 11 and are detuned from the ^^Rb D2 line by about 



IX. SUMMARY 

This paper discusses spectral and transport properties of the 
cold atom analog of a 2DEG in a lattice, subject to a non- 
Abelian gauge field with C/(2) symmetry. 

In the continuum limit of the lattice, the system maps onto 
two oppositely-charged coupled harmonic oscillators, with 
a coupling constant proportional to the strength of the non- 
Abelian field. The Landau energy levels of the Abelian prob- 
lem evolve into entangled states of this particle-hole pair. 

These features also characterize the energy spectrum of 
the the corresponding lattice problem. In fact, the transi- 
tion from Landau levels to Landau bands is the analog of the 
generalization from the butterfly to the moth spectrum as the 
Abelian system becomes non-Abelian. The non-Abelian cou- 
pling breaks the degeneracy of the Landau levels; the spec- 
trum depends explicitly on the transverse momentum. 

The non-Abelian system exhibits antiferrimagnetic-type 
ground states, whose components, A and B, need not be de- 
generate, and in fact may have very different transport prop- 
erties. A particularly interesting example of this is the zero- 
energy state for a2 = 1/4, where the degenerate A and B 
components have different localization properties. Addition- 
ally, an intriguing relationship between the A and the B com- 
ponents occurs for 02 = ^, as these two components cor- 
respond to the positive and the negative energy states of the 
system. Such novelties may open new avenues for exploring 
frontiers of physics with cold atoms. 
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The use of ultracold atoms to simulate relativistic as well 
as non-relativistic theories and study the effect of disorder is 
an exciting field of research. In a two-dimensional lattice sub- 
ject to a non-Abelian gauge field, one can induce not only lo- 
calization transitions, but also a transition from relativistic to 
non-relativistic theory by tuning the lattice anisotropy. A well 
known feature of the Dirac Hamiltonian is an extra term in 
the conductivity attributed to Zitterbewegung (ZB) II16I1 cor- 
responding to inter-band transitions. It has been suggested 
that such a term is responsible for the finite conductivity of 
graphene described by a massless Dirac energy spectrum lllql . 
In other words, it is ZB that makes it impossible to localize rel- 
ativistic particles, as it is connected with the uncertainty of the 
position of a relativistic quantum particle due to the creation 
of particle-antiparticle pairs. Therefore, the origin of der- 
ealization characterizing the non-Abelian system that persists 
even for infinite disorder (A — > 06) can be attributed to ZB. 

The detection of relativistic particle and a transition from 
non-relativistic to relativistic dispersion in cold atoms in opti- 
cal lattices was recently discussed; it was shown that the rel- 



ativistic dispersion can be detected using atomic density pro- 
files as well as Bragg spectroscopy 11 3i1 . 

Our detailed study for various values of 02 captures some 
of the universal features of non-Abelian systems. Exploration 
of the two-dimensional space (ai, a2) may reveal additional 
phenomena, and the richness of U{N) gauge systems with 
N > 2 remains to be explored. Moreover, the effects of inter- 
particle interactions remain to be investigated [|19|1 . 
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